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ABSTRACT 

We present a detailed analysis of the XMM-Newton and Chandra X-ray data of the young Type la supernova remnant SNR 0519- 
69.0, which is situated in the Large Magellanic Cloud. We used data from both the Chandra ACIS and XMM-Newton EPIC MOS 
instruments, and high resolution X-ray spectra obtained with the XMM-Newton reflection grating spectrometer (RGS). 
Our analysis of the spatial distribution of X-ray line emission using the Chandra data shows that there is a radial stratification of 
oxygen, intermediate mass elements (IME) and iron, with the emission from more massive elements peaking more toward the center. 
Using a deprojection technique we measure a forward shock radius of 4.0 ± 0.3 pc and a reverse shock radius of 2.7 ± 0.4 pc. 
We took the observed stratification of the shocked ejecta into account in the modeling of the X-ray spectra, for which we used multi- 
component non-equilibrium ionization models, with the components corresponding to layers dominated by one or two elements. An 
additional component was added in order to represent the shocked interstellar medium, which mostly contributed to the continuum 
emission. This multicomponent model fits the data adequately, and was also employed to characterize the spectra of distinct regions 
extracted from the Chandra data. From our spectral analysis we find that the approximate fractional masses of shocked ejecta for the 
most abundant elements are: Mo « 32%, Msi/s ~ 7%/5%, MA,+ca ~ 1% and Mp e » 55%. From the continuum component we derive 
a circumstellar density of n H = 2.4 ± 0.2 cm 4 . This density, together with the measurements of the forward and reverse shock radii 
suggest an age of 0519-69.0 of 450 ± 200 yr, somewhat lower than, but consistent with the age estimate based on the extent of the 
light echo (600 ± 200 yr). 

Finally, from the high resolution RGS spectra we measured a Doppler broadening of a = 1873 ± 50 kms~', from which we derive 
a forward shock velocity of v FS = 2770 ± 500 kms~'. We discuss our results in the context of single degenerate explosion models, 
using semi-analytical and numerical modeling, and compare the characteristics of 0519-69.0 with those of other Type la supernova 
remnants. 

Key words. X-rays: individuals: SNR05 19-69.0 — ISM: individuals objects: SNR05 19-69.0 — ISM: supernova remnants — method: 
data analysis 



1. Introduction 

Thermonuclear (type la) supernova explosions have drawn a lot 
of attention over the last decade, as they provide a powerful 
tool to measure cosmological distances (Perlm utter et al.fl999| 
Ries s et aL]|1998[ ). The reason is that they are bright enough 
to be observed over large distances, and that, compared to core 
collapse supernovae, their peak luminosity shows relatively lit- 
tle variation, which can be further reduced by applying the 
empirical peak luminosity — decline rate correlation (Phillips 
et al.|1999]>. Ne vertheless, some evolutionary systematic effects 



(Panagia 2005) may take place, which may affect the measured 



values of cosmological parameters. Thus, it is important to un- 
derstand the mechanism and physics that governs thermonuclear 
explosions. 

In addition to the extensive studies of the supernovae them- 
selves, through their light curves (e.g |Woosley et al.| 2007 



reverse shocks. The reverse shock, propagating inwards into the 
supernova ejecta, efficiently heats the metal-rich matter. The hot 
plasma produces X-ray spectra, abound with prominent emission 
lines. These spectra carry the imprints of the chemical composi- 
tion of the plasma, and, therefore, the distribution of the elements 
in the supernova ejecta, which are determined by the explosion 
properties. Thus, a detailed analysis of the X-ray data of super- 
nova type la remnants provides a powerful complementary tool 
in the studies of the nature of a thermonuclear supernova pro- 
genitor (e.g. Sorokina et al.|2004 Badenes et al.|2003) . 

The current generation of X-ray observatories, such as 
XMM-Newton and Chandra, provide the capability to obtain si- 
multaneously spectral and imaging data of extended objects such 
as SNRs. Young SNRs in the Large Magellanic Cloud (LMC) 
are especially suitable targets for these telescopes, as the rel- 
ative proximity (48 kpc, e.g. Pietrzynski et al. |2009 Koerwer 



Blinnikov et al. 2006) and spectra (e.g. |Branch et al.|2009[ ), one 
can also turn the attention to the SN type la remnants. A typical 
young supernova remnant (SNR) is a bright X-ray source, due to 
the high temperatures of the plasma, heated by the forward and 



2009 ) of the LMC yields a sufficiently high signal-to-noise level 
of the X-ray data as observed by the CCDs and grating spectrom- 
eters on board Chandra and XMM-Newton. The SNRs are still 
large enough (typically ~ 1') to perform morphological stud- 
ies with Chandra and small enough to be excellent targets for 
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ods and techniques we used to deal with the spectral and imag- 
ing data are presented in Sect. [3] Sect. [4] contains the principal 
results of the study, which are discussed in Sect. [5] We conclude 
the paper in Sect. [6] 



2. The data overview 



Fig.l. Chandra smoothed RGB image of SNR 0519-69.0. Red 
— 0.5 - 1.0 keV, green — 1.0 - 2.0 keV, blue — 2.0 - 4.0 keV. 
North is up, east is left. 



the Reflection Grating Spectrometer of XMM-Newton. An ad- 
ditional advantage is the relatively low interstellar absorption 
column towards LMC SNRs {Nh ~ 10 21 cm -2 ), as compared 
to most Galactic SNRs. 

In this paper we present an X-ray study of the LMC 
SNR 0519-69.0. The SNR has an irregular patchy morphology 
and an angular extent in X-rays of 33", corresponding to a shell 
radius of 4 pc. A three color X-ray image of SNR 0519-69.0, 
based on the Chandra data, is presented in Fig. [T] 



In the optical SNR 0519-69.0 has been investigated by Tuohy 
|etaLl l fl982) , [Smith et al.| ( fl99TT ) and |Ghavamian et al.| l |2007p t 
is also one of several LMC SNRs for which a light-echo has been 
identified, from which an age of 600 + 200 years can be deduced 

( [Rest et al.|2005) . 

In X-rays the remnant was studied by Hughes et al. ( 1995[ l, 
who analyzed the ASCA spectra and found that the SNR is 
oxygen-poor and iron-rich and must be a remnant of a thermonu- 
clear supernova (SN la) explosion. The analysis of Chandra data 
by | Williams et al.| (|2001|l revealed the separation between the 



shocked ejecta and the shocked circumstellar medium (CSM). 

In the present study we analyzed archival X-ray data of 
SNR 0519-69.0 from the both XMM-Newton and Chandra ob- 
servatories. We used the ACIS data of Chandra and the imag- 
ing spectroscopy (EPIC) and the high spectral resolution grating 
(RGS) data of XMM-Newton. 

Various techniques to analyze the available data are em- 
ployed. The high spatial resolution Chandra data provide us with 
an opportunity to study the composition of the SNR as a func- 
tion of radius. Whereas the RGS data offer high resolution spec- 
troscopy, but without much spatial information. 

For the analysis of the SNR 0519-69.0 X-ray spectra, we 
employed the spex (Kaa stra et al.|1 996) spectral fitting software 
(version 2.01.05, November 16, 2009). The package contains the 
most up-to-date atomic data and has a wide range of plasma 
emission models, which is especially helpful and important in 
fitting complex spectra from objects such as SNRs. 

We fit the XMM-Newton EPIC and RGS spectra with single 
and multicomponent NEI models. The high spectral resolution 
RGS data allow us to resolve details of the Fe-L line emission, 
and to measure line velocity broadening, due to the thermal and 
bulk motion of the shocked supernova ejecta. 

This paper is organized as follows. First in Sect. [2] we de- 
scribe briefly the Chandra and XMM-Newton data. The meth- 



2.1. Chandra 

SNR 0519-69.0 was observed only with the ACIS instrument of 
Chandra X-ray Observatory (obs ID 1 18) for 41.1 ks on June 21, 
2000. We analyzed the data with CIAO 4.1.1 software product 
and CalDB 4. 1 .2 calibration data. Apart from an X-ray spectrum 
of the whole SNR, we also used our own software to extract 
spectra from certain regions based on emission characteristics. 
The software reads the standard pipeline event lists and selects 
events from it based on a "mask" image, i.e. an image containing 
only zeroes and ones. An event is selected for spectral extraction, 
if its sky coordinates corresponds with a pixel of value one in the 
mask image. The event grade selection and spectral data binning 
scheme are identical to the standard CIAO software. The advan- 
tage of this method is that spectra can be extracted from regions 
with more complex shapes than with the region files, and can 
be easily constructed based on hardness ratios or narrow band 
images. These spectra were used together with the ancillary and 
main response files generated for the SNR as a whole. Given that 
the remnant occupies only a small part of the ACIS-S3 CCD, us- 
ing these files for the response is justified. For the background 
spectrum, we extracted a spectrum from an annulus with radii 
spanning from 19" to 37", centered on the SNR. 



2.2. XMM-Newton 

We used XMM-Newton observations of SNR 0519-69.0 (obs 
ID 0113000501) that were made on September 17, 2001, and 
have a total exposure of 47.8 ks. 

For our study we concentrated on the EPIC MOS and RGS 
data. Although the EPIC MOS instruments ( |Turner et al.|2001| l 
have a lower sensitivity than the EPIC pn instrument, they have 
a higher spectral resolution, which is important for line-rich 
sources, such as SNRs. The MOS1 observations were performed 
with thick (25.4 ks) and medium (15.6 ks) filters, whereas the 
MOS2 data were acquired with the medium filter only. For the 
spectral analysis the spectra, response and ancillary response 
files were combined, using weights proportional to the exposure 
times. For the background, spectra were extracted from an annu- 
lus with an inner radius of 2.1', and an outer radius of 4.3'. This 
is larger than for the Chandra background spectrum, because the 
XMM-Newton point spread functions have much broader scat- 
tering wings. 

The RGS is a slitless spectrometer ( den Herder et al.|2001| l. 
For an extended source, this implies that the spectrum is smeared 
by the image of the source itself. For Small/Large Magellanic 
Cloud remnants, the smearing is modest, but present, and it gives 
rise to a change in the line spread function. For our analysis, 
we incorporated this effect into the response matrix by convolv- 
ing the standard (point source) response matrix with the bright- 
ness profile of the SNR, as obtained from the Chandra obser- 
vations. All the results presented here made use of these modi- 
fied matrices. This procedure was also applied to the RGS data 
of the SN 1006 remn ant ( |Vink et al.|2003] > and SNR 0509-67.5 
( |Kosenko et al.|2008) . 
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Apart from adapting the RGS response matrix, all reduc- 
tion for both MOS and RGS data were made with the standard 
XMM-Newton software package SAS version 7.1.0. 



3. Data analysis 

3. 1 . Chandra radial emissivity profiles 

The high spatial resolution Chandra data allow us to plot radial 
emissivity profiles of the remnant's shell, with a deprojection 

( [200% For 
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technique that was also employed by|Helder & Vink| 



the measurements, we first made a radial surface brightness pro- 
file of the remnant, using the centroid of the broadband emis- 
sion (0.5-8.0 keV) as center. Then, assuming that the SNR is 
spherically symmetric, we used the Lucy-Richardson technique 
( Lucy|[l974} |Richardson fl972| ) to deproject the surface bright- 
ness profile into an emissivity profile (following the method of 
Willingale et~aT1[T996] |Helder & Vink|[2008| . We applied this 
technique to images in different energy bands, where the emis- 
sion lines from different species dominate. For the oxygen band, 
we choose 0.5-0.7 keV, for the added silicon and sulfur band, we 
took 1.75-3.00 keV. Because the spectra indicated that in the Fe- 
L band (~ 0.7 -1.1 keV) the relative contribitions from different 
Fe-L shell ions varies as a function of position (Sect. 3.5.2| i, we 
divided the iron energy band in a low (0.7-0.9 keV) and high 
(0.9-1.1 keV) energy band, which correspond roughly to line 
emission dominated by Fe XVII and Fe XX ions respectively. 

For the broad band (0.5-8.0 keV) energy range, we calcu- 
lated an inner and outer radius. The inner radius was determined 
where the emissivity (i.e. deprojected, e) drops below one half 
of the maximum value. We ignored emissivities at radii less 
than 7", since the deprojection technique is not well constrained 
at these lower radii, due to the low contribution to the surface 
brightness X of shells at small radii (E oc e(R)R 2 ). In general, 
there is no sharp rise in emissivity at the outer radius, so taking a 
value of one half of the maximum would underestimate the outer 
radius. Hence, we took a value of 1/20 of the peak value, which, 
in the deprojections, is just above the noise level of the back- 
ground at large radii. In the remainder of this paper we assume 
that this measured inner radius of 2.7 ±0.4 pc, corresponds to the 
location of the reverse shock and the outer radius to the forward 
shock at 4.0 + 0.3 pc. We estimated the error from the standard 
deviation of this measurement repeated for 18 individual slices 
of 20° each. 

We use the deprojected profiles to determine the shells in the 
remnant in which different species dominate (Fig. [2]), which we 
will use in section [4] for estimating total masses of the shocked 
species. Table [T] lists the ranges in which the different compo- 
nents dominate. Silicon and sulfur dominate the emissivity in 
region 2/3, however, our numerical models (see Sect.{5]l show 
that the Si/S and Ar/Ca layers coincide. For this reason, we split 
region 2/3 evenly in two parts. Furthermore, we chose 3.6 pc as 
the division between oxygen and the shocked CSM, since for 
this value, the volume of the shocked CSM is one fourth of the 
volume of the total remnant. The latter seems reasonable, if we 
assume a compression ratio of 4 over the shock front during the 
whole lifetime of the remnant. Note that all values mentioned 
above contain uncertainties. For the corresponding volumes, we 
estimate errors of 50%. In Fig. [2] we also show the observed ra- 
dial profile and its deprojection for the 5-7 keV hard X-ray band. 
This band includes both continuum emission and Fe-K shell line 
emission around 6.5 keV. Unfortunately, the signal to noise ratio 
in this band is poor, but the profile suggests tentatively that the 
5-7 keV emission is associated mostly with the Fe XX line emis- 
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Fig. 3. Combined fitting of the EPIC MOS (scaled down by a 
factor of 10) and ACIS (scaled up by a factor of 10) data with 
single ionization timescale NEI model of spex. 



sion in the 0.9-1. 1 keV band. With some caution, it also suggests 
that there is no strong non-thermal X-ray continuum associated 
with the forward shock. 



3.2. EPIC and ACIS spectra 

Fitting the XMM-Newton EPIC MOS and Chandra ACIS spec- 
tra with single-ionization timescale non-equillibrium ionization 
(NEI) spex model provides us with the typical values for the 
emission measure, electron temperature, ionization timescale 
and abundances in the remnant. The corresponding spectra and 
the best-fit model are shown in Fig. [3] T he best-fit parameters of 
the spectral model are listed in Table |2J The best-fit abundances 
are plotted in Fig. |4] For comparison, also the abundances, de- 
rived from the analysis of the younger supernova Type la rem- 
nant 0509-67.5 ( jKosenko et al.|2008| > are plotted. 

A single ionization timescale NEI model is not the best ap- 
proach to explain emission from such a complicated and lay- 
ered object as a SNR, because in the shell one expects ionization 
timescale and abundance gradients. For example, Fig. [3] clearly 
shows that Si Lya line is missing in the model and Fe K line flux 
is overestimated. This discrepancy leads to a very high value of 
the fit statistic; x 2 /d.o.f. ^ 30. 

Apart from the MOS, we fitted also the EPIC pn data with 
an NEI model, that yields approximately the same values of the 
basic parameters of the spectrum. 



3.3. XMM-Newton RGS spectra 

The RGS spectra of SNR 0519-69.0 were also fitted with spex 
NEI models. The fit range was limited to the 0.5-1.1 keV range, 
which has the more prominent line emission. The best-fit param- 
eters of the single NEI model are listed in Table [2] (first column). 
This single-ionization timescale model was not able to repro- 
duce the fluxes of O VII and Fe XXI ions: the modeled lines 
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Table 1. The inner and outer radii for the SNR 0519-69.0 shell in different energy ranges. 
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0.5- 


8.0 keV 


CSM 


O 


Si/S 


Ar/Ca 


Fe-high 


Fe-low 











1 


2 


3 


4 


5 


Rm (pc) 


2.7 


+ 0.4 


3.6 


3.48 


3.405 


3.33 


3.15 


2.7 


^out (PC) 


4.0 


+ 0.3 


4.0 


3.6 


3.48 


3.405 


3.33 


3.15 



X! 
CD 
N 

s 

u 
o 
c 



R [arcsec] 
12 14 1 



R [arcsec] 
10 12 14 16 



1.0 
0.5 
0.0 
-0.5 




0.50-0.70 keV 

0.70-0.90 keV 

0.90-1.10 keV 

1.75-3.00 keV 

5.00-7.00 keV 



-0.5 



-a 

D 
N 

6 
o 

IJJ 



2.0 2.5 



3.0 3.5 
R [pc] 



4.0 2.0 2.5 



3.0 3.5 
R [pc] 



4.0 



Fig. 2. Azimuthally averaged radial profiles in different energy bands. Left panel: surface brightness radial profiles. Right panel: 
deprojected emissivity profiles. The vertical light grey lines correspond to the spectral components, used in Sect. 3.4 



Table 2. The best-fit NEI parameters of the XMM-Newton EPIC, RGS and Chandra ACIS data. Columns 1 and 2 list the data for 
the entire remnant, column 3 lists data for the inner shell of the SNR and column 4 — for the outer shell (see Fig.|7]i. 



Parameter 


RGS 


EPIC MOS + ACIS 


ACIS inner 


ACIS outer 


n e n„ V (10 58 , cm" 3 ) 
kT e (keV) 
n e t (10 10 , s cm" 3 ) 
n„ (10 21 , cm" 2 ) 


i-8-fi 

2 2i+ - 08 
2 69"^ 


q 7+0.2 
J - -0.2 
9 Q9+0.03 

2 27+0 01 
2 62+^ 


2 63 +0 33 

1 CQ+0.06 

' -R-AS) 
9 99+0.02 

-Q.0Q2 
9 74+0.11 


66 +6 45 

u,uu -0.08 

1 7^+0.14 
1,7 -1.07 

3.90 ± AT/A 


X'ld.o.f. 


2.4 


29.4 


10.4 


2.5 



are weaker than the data indicate (top panel of Fig. [5]). For this 
reason, and inspired by the spatial layering, as indicated by the 
Chandra data (Sect. 



3.1 



we fitted the RGS spectra in the range 
of 0.5-1.1 keV with a three-component NEI model: one compo- 
nent is for pure oxygen, the second — for the low-ionized iron 
(Fe XVII - Fe XVIII), and the third — for the high-ionized iron 
(Fe XIX - Fe XXI). The corresponding spectra are shown in the 
bottom panel of Fig. [5] The corresponding best-fit parameters 
are listed in the Tablej3] As can be seen in Fig. [5] this model 
gives a much better fit to the data than the single NEI model. 
Note, that no Ne emission is needed to obtain a good fit (Ne IX 
has a prominent line at 0.92 keV and Ne X at 1.02 keV). 



Table 3. The best-fit parameters of the RGS spectral fitting with 
the three-component NEI model. x 2 /d.o.f. = 1.8 



Parameter 



O 



Fe-low Fe-high 



n e n„ V (10 58 , cm" 3 ) 
kT e (keV) 
n e t (10 10 , s cnr 3 ) 
orgs (km s" 1 ) 



50+ 08 

9 ^+0.8 
^'"'-0.5 



5.6-8;| 4.8!i; 2 

i 99+O.O5 9 c i +6.54 

1,z -0.33 z " JL -0.99 

3 +23 5 2 +34 

J - -0.2 J - -1.4 

1873 ± 50 



An important advantage of the high spectral resolution RGS 
spectra is that it enables us to measure the line broadening. The 
measured best-fit value of the line velocity broadening is cr v = 
1680 + 50 km s" 1 for the single NEI model and cr v = 1873 ± 50 
kms -1 for the multi-component model. We adopt the last value 



as the most reliable value for the Doppler broadening, since the 
model fits the data much better (Fig. [5]). 
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Fig. 4. Best-fit abundances of SNR 0519-69.0 and SNR 0509- 
67.5 (Kosenko et al. 2008[ l, derived from single-ionization 
timescale NEI models. The data are in solar units (An ders &| 
Grevesse|1989| l, normalized by Si abundances. 
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3.4. EPIC MOS and RGS spectra combined, 
multicomponent approach 

A next step in our study was to combine both the EPIC MOS 
and RGS data and try to model it with a multi-component NEI 
model, in which each component roughly corresponds to the lay- 
ering observed in the radial profiles (Sect. [XT]). In addition to the 
three components used for fitting the RGS spectra (Sect. 3.3 1, 



we added three additional components, which contribute mostly 
to the emission outside the RGS spectral range. In total this 
amounts to six components: oxygen (O, or component 1), sil- 
icon/sulfur (Si/S, or component 2), argon and calcium (Ar/Ca, 
or component 3), high ionized iron (Fe-high, or component 4), 
low ionized iron (Fe-low, or component 5), and a continuum 
dominated model (component 0), for which we used an NEI 
model with the abundances fixed to those of the LMC ( |Russell| 
|& Do pita 1992). In addition the model included an interstellar 
absorption component. 

The corresponding principal abundance parameters in com- 
ponents 1-5 were set to 10 7 times the solar value for the spe- 
cific elements considered, whereas the abundances of other el- 
ements were set to zero. Thereby we secure that the absolute 
abundances of the corresponding elements are at least two or- 
ders of magnitude higher than those of the hydrogen, i.e. for all 
practical purposes these components correspond to pure metal 
plasma's. The corresponding spectra and the best-fit multicom- 
ponent model are presented in Fig. [6] The best-fit parameters 
(emission measure, electron temperature, ionization timescale) 
are listed in Table [4] The best-fit value for the hydrogen column 
density is 1.7^ x 10 21 cirT 2 . 

The multi-component NEI model gives an adequate fit to 
the data. In contrast to the three NEI component model used 
to fit only the RGS data, some of the oxygen line emission in 
this multi-component model attributes part of the corresponding 
line emission to the shocked CSM component, in particular the 
O VIII line emission. One peculiarity of the multi-component 
model is the high n e t value for Ar/Ca. This is not an artefact 
of the model chosen. Isolating that part of the spectrum that is 
dominated by Ar/Ca line emission and searching for other com- 
binations of kT e and n e t resulted in similar best fit values. In 
general, « e f correlates inversely with kT e . This means that if the 




0.8 

Energy (keV) 

Fig. 5. Top panel: the four RGS spectra RGS1, RGS2 of 
SNR 0519-69.0. Shown are both first order (upper frame) and 
second order (lower frame) spectra and best fit single NEI model 
including a line velocity broadening model. Bottom panel: The 
same as the top panel, but the spectral model consists of three 
NEI components. Black crosses are the data, solid lines are the 
models. 

n e t of the plasma is in reality lower than our best fit value, an 
unrealistically high electron temperature is required. 

3.5. Spatially resolved spectroscopy 

3.5.1 . Inner and outer rings of the remnant 

In the analysis of the Chandra data we extracted ACIS spectra 
from two shells: the outer one which is presumably the shocked 
CSM of the remnant and the inner one which is the region where 
the emission is produced by the shocked ejecta material. These 
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Table 4. Parameters of the components in the multi NEI approach. Errors are lcr rms, Ax 2 = 2,x 2 Id.o. f = 7.8. Emission measure 
EMx - n e rixV. 

Component number 1 2 3 4 5 

Layer O Si S Ar Ca Fe-high Fe-low CSM 



(10 54 cm" 3 ) 8.87^ 2.47!"' 1.56^ 0.78!°,« 0.60^ 2.52*}fj 2.15$" 2.36$» X 10 5 

*T„(keV) 0.72 + |j : ™ 7.00+™ 7.00+™ 2.67+1™ 2.67 + |™ 2.73 + ™ 1.26 + °f 0.64^ 

M(10 10 scm3) 1.52?** 3.63*$ 3.63^ 27$™ 27.2 + _± 3.78*& 2.55+° : « 80.0% 

n e /n x 6.6 21.2 33.7 30.5 39.6 19.6 17.5 1.2 



0.001 




0.6 0.7 0.8 0.9 1 

Energy (keV) 

Fig. 6. The EPIC MOS 1 and MOS2 combined (upper frame) and 
the RGS 1 and RGS2 first order (lower frame) spectra with the 
best-fit model of the six NEI components and interstellar absorp- 
tion. Black crosses are the data, the solid lines show the model. 
The contributions from the different components are marked 
with the numbers listed in Table [4] The dotted lines correspond 
to the CSM component (0). 



regions are outlined with green circles in Fig. [7] Each of the 
spectra was analyzed in the usual way and fitted with single 
NEI models. The corresponding best-fit parameters are listed in 
Tabled 



3.5.2. ACIS oxygen and iron masks 

The low energy band (0.5-1.1 keV) of the ACIS spectra was 
treated more thoroughly. We extracted three images of the rem- 
nant in the energy bands (see the RGS spectra with the corre- 
sponding emission lines) of 0.5 - 0.7 keV (O only), 0.7 - 0.9 
keV (Fe-low) and 0.9 -1.1 keV (Fe-high). Next, three spatial 
masks were created on the basis of these images, in order to lo- 
cate regions with the relative emission peaking in either of these 
energy bands. The masks are presented in the Fig. [7] They were 
applied to the ACIS image of the remnant. Finally, we extracted 
spectra using these, non-overlapping, image masks. 




Fig. 7. Left panel: Chandra image of SNR 05 19-69.0. Outer and 
inner regions that were considered are outlined with green lines. 
Right panel: the masks of the Chandra image of SNR 0519-69.0: 
red — oxygen (0.5 - 0.7 keV), green — Fe-low (0.7 - 0.9 keV), 
blue — Fe-high (0.9 - 1.1 keV) 



The resulting spectra are poorly fitted with single-ionization 
timescale NEI models, the values of the best-fit abundances con- 
tain large errors, and physical parameters, such as temperature 
and ionization timescale are of the order of the typical values 
obtained in the analysis of the EPIC and ACIS data of the entire 
remnant. This suggests that the masks did not totally separate out 
specific layers, but rather that certain layers may dominate the 
emission in certain regions, with some overlap from other lay- 
ers. Therefore, we applied the six-component NEI model with 
the parameters obtained in the fitting of EPIC and RGS data 
to each of the masked spectra. We fit only emission measures 
(i.e. normalizations) of each component in order to evaluate the 
contributions of the different components, each corresponding 
to certain elements, to the spectra of these three distinct regions. 
The spectra and the best-fit models are presented in Fig. [8] 

The contributions of each component to the different spec- 
tra are presented in the Fig. [9] in which we plot the emission 
measures of each NEI component. 

4. Results 

The fitting of the XMM-Newton EPIC MOS and Chandra ACIS 
spectra with spex NEI models allowed us to measure (Table [2} 
the parameters of the plasma, such as typical electron tempera- 
ture kT e a 3 keV, ionization timescale n e t = 2.3xl0 10 s/cm and 
abundances (plotted on Fig. [4]) in the remnant. The derived best- 
fit abundances are similar to those of SNR 0509-67.5 (Kosenko 
|et al.||2008) , but overall the 0519-69.0 remnant is more metal- 
rich (S, Ar, Ca and in particular Fe) as it is older and more of the 
supernova ejecta matter has been swept up by the reverse shock. 
The relative amount of the light elements (such as O, Ne, Mg) is 
comparable for both SNRs. 
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4.2. Structure of the ejecta 

The radial emissivity profiles of the SNR in different energy 
bands, drawn from the Chandra observations, clearly show strat- 
ification of the elements in the remnant. Fig. [2] shows that the 
outermost layer is oxygen rich, the next inward layer produces 
most of the Si and S emission, all the iron emission comes from 
the innermost layers of the shell. The iron shell is split in two 
regions with different ionization properties. The outer layer of 
high ionized (Fe XX) and the inner one of low ionized (Fe XVII) 
iron reflect the effects of time-dependent ionization processes in 
the shocked ejecta: the inner layer has been shocked later by the 
reverse shock, and therefore has the lowest n e f . 

This stratification is also confirmed by the analysis of the 
masked Chandra images (Fig. [7] right panel) and their spectra 
(Fig. [§). Fig. [9] shows the contribution of different elements to 
the emission from the masked images. 



Fig. 8. The three masked ACIS spectra. From top to bottom: 
oxygen (data shifted upwards for three orders), Fe-high, Fe-low 
(data shifted downwards for three orders). 
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Fig. 9. The emission measure of each of the six NEI components 
shows the contribution of different elements to the emission of 
the masked regions of the Chandra data. 



4.1. The circumstellar matter density 

The multicomponent NEI model provides us with the parame- 
ters of the CSM component (component in TableQ. Using the 
best-fit emission measure value we can make an estimate of the 
density of the CSM in the vicinity of the remnant. 



1 



«CSM - 



EM; 



LMC 



2 V VsNR«f/«H 



= 2.4 ± 0.2 ctrT 



(1) 



where £M LM c = 2.36+" ?" x 10 59 cm 3 and n e /n H = 1.2 are 



-o.n 



taken from Table [4] We assumed, that the volume of the shocked 
CSM is 1/4 V SNR , with V SNR = 47T/3 R 3 out = 8.0 x 10 57 cm 3 and 
R om =4.0 pc (Table [T}. 



4.3. The chemical composition 

The measurements of the inner and outer radii of the shell in 
different energy bands (Table [TJ allow us to estimate the masses 
of the layers, which have been swept up by the reverse shock. 
Using the best-fit emission measure value EM X = n e rixVx of 
the each model of the multicomponent fitting (Table [4}, we can 
express the mass M x of an element "X" (i.e. O, Si, S, Ar, Ca, 
Fe) as 



M x = n x V x m x 



[em 

V * 



EM X V X 
/«x 



OT U A > 



(2) 



where n e /n x — electron to ion densities ratio, m v — atomic 
mass unit, A x — atomic mass of the element, V x — volume 
occupied by the element (derived from the radii of Table [TJ. 

Combining the available data we obtain Mo = 
O.36M , Msi/s = 0.14/0.10M o , M Ar/Ca 
0.08/0.07 M Q , M Fe = 1.05 M Q , M tot = 1.8 M Q . Note, that 
the mass values for the elements are affected by the ~ 50% 
errors of the volume estimates, although the total mass is less 
affected. These masses are higher than may be expected for a 
single degenerate Type la SN model, for which the total mass 
should be M tot = 1 .4 M . The most likely reason for this discrep- 
ancy may be the errors in the volume estimates or in the filling 
factor. In particular, inhomogeneities such as knots may give 
emission measures skewed toward the higher density regions. 
In fact, the supernova material may be porous and clumpy, thus 
the emitting volume should be scaled down using a filling factor 
(e.g. see 3D simulations of thermonuclear explosions by Ropke 
|& Bruckschen|2008] l. For example, if we assume volume filling 
factor in the supernova of ~ 0.4, then we get more adequate 
values M = 0.23 M Q , M si/S = 0.09/0.07 M Q , M M/Ca = 
0.05/0.05 M , M Fe = 0.67 M , M tot = 1.15 M . 

The multicomponent fitting of the XMM-Newton EPIC and 
RGS spectra reveals a specific behavior of the different NEI 
components. Table |4] shows that emission measure (EM X = 
n e tixVx) of each component drops from oxygen to calcium, 
while the values of the ionization timescale parameter n e t rise 
from O to Ca. This trend indicates that increasing from the 
lighter to the heavier elements, the electron density overcomes 
the factor of time, which probably does not change drastically 
for the shocked intermediate mass elements (from Si to Ca). 
Thus, qualitatively, the decrease of the emission measure and in- 
crease of the ionization timescale reflect the relative abundances 
of the species from O to Ca, i.e. the quantity n e rixVx/(n e t) = 
nxVx/t cc Mx (total shocked mass of an element X) is lower 
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for the heavier species. A similar behavior of the ionization 
timescale parameter for different species was established by 
Lewis et al. (2003) in their analysis of the Chandra data of 
SNRN103B. 

From these considerations, applying the time of t 500 
years to all components, we can make rough estimates of the 
relative contributions of the species in the SNR 0519-69.0 ejecta 
to be as follows: M = 32%, M si/S = 7%/5%,M Ar+Ca = 
1%, Mp e = 55% of the total shocked mass, which are roughly 
consistent with the mass ratios derived above, using volume es- 
timates. The relative amount of Ar and Ca is probably underes- 
timated, due to the high value of ionization timescale parameter 
measured in the component 3, for which we do not have a good 
explanation (note the corresponding errors bars in Table|4|. 



5. Discussion 

5.1. Stratification and composition of the ejecta 

Using the Chandra images we built azimuthally averaged ra- 
dial emissivity profiles of the SNR in different energy bands, 
which clearly demonstrate the layering of the elements in the 
supernova (Pig. |2j. Also these profiles reveal the effects of time- 
dependent ionization, as the emission from low ionized iron 
originates mostly from the inner layers of the iron shell, which 
is distinct from the high ionized "older" outer iron layers. 

These effects are also visible in the Chandra spectra extracted 
using the masked images (Fig. [7J. However, apparently projec- 
tion effects play a role in that the spectra of the more inward ly- 
ing layers are contaminated by the outer layers. For example, the 
spectral fitting with multi-component NEI models shows (Fig.|9| 
that the oxygen-rich component (component 1) is equally abun- 
dant in all three (O, Fe-high, Fe-low) masks. Si and S (compo- 
nent 2) is almost absent in the O-mask, as it is inside the oxygen- 
rich layer, but it does contribute to the spectra of inner layers 
of the Fe-low mask. Ar and Ca (component 3) tend to be in the 
outer layers, but the parameters of this component cannot be reli- 
ably measured, and have considerable errors. The O-mask spec- 
trum hardly contains any contribution from the iron-rich com- 
ponents (4,5), as this outer layer can be isolated without much 
projection effects. 

Fe-K emission is present in both the regions corresponding 
to the Fe-low and Fe-high masks. This is probably due to pro- 
jection effects, since our spectral analysis suggests that the Fe-K 
emission is mostly coming from the high ionization layer, which 
was shocked earlier (component 4, Table [4j. The radial emis- 
sion profiles are consistent with this association, but we note 
that the signal to noise of the 5-7 keV image on which this is 
based is low. The high temperature of this high ionization layer 
may be surprising, but it should be mentioned that the electron 
temperature depends on both the shock velocity, and the degree 
of equilibration between electron and ion temperature. In gen- 
eral, one expects that ejecta that is shocked later by the reverse 
shock is heated to higher temperatures, but this may be off set 
by the shorter time available for the electrons to equilibrate with 
the ions. The situation in 0519-69.0 may somewhat resemble the 
situation in the south of Kepler's SNR. The analysis of |Cassan>] 
Chenai' et al. ([2004) shows that for this region Fe-K peaks at a 



5.2. Comparison of the stratification with other Type la 
supernova remnants 

The chemical stratification of SNR 0519-69.0 is striking. Even 
more so if one considers that this SNR is much more distant than 
the well studied Galactic Type la SNRs Tycho (SN 1572), Kepler 
(SN 1604) and SN 1006. All these SNRs also show some strat- 
ification, but less pronounced. In some cases the Fe, Si, and O 
layers seem to be very close together or even indistinguisable as 
separate layers. The latter is especially true for Tycho, in which 



silicon and iron are found very close to the shock front ( Warren 
et al.|2005] Badenes et al.|2006 1. For Kepler the situation is dif- 



ferent for the northern part of the SNR than for the southern part 
( |Cassam-Chenai et al.|20 04). In the north the spatial separation 
between Si and Fe is small, whereas in the south the situation is 
in fact similar to the stratification of SNR 0519-69.0. 

The situation for SN 1006 is a bit more complicated. There 



the O-rich ejecta lie very close to the forward shock (Cassam- 
|Chenai et al.|2008| ), but otherwise the Fe-rich layer seems to lie 
at smaller radius, and may not even have been shocked by the 
reverse shock, since UV absorption spectra shows the reverse 
shock to go through Si-rich material (Hamilton e t al.|2007] ). 

So what is the origin of this variation in stratification among 
Type la SNRs? It could be that there is some variation among 
Type la supernovae itself, with some supernovae producing more 
mixing than others. Optical spectroscopy of Type la supernovae 
suggests a strong radial stratification (e.g. Bra nch et al.| |2005 
Mazzali et al. 2007 ) with the the outer ejecta consisting of low 
mass elements (O, Mg), then followed by intermediate mass el- 
ements (IME), and the inner ejecta consisting of Fe and Ni iso- 
topes. However, a group called High Velocity Gradient (HVG) 



larger radius than the Fe-L emitting region. 



supernovae ( Benetti et al.||2005) show a strong gradient in the 
Si II velocities as a function of time, and evidence for high ve- 
locity I ME. For the HVG S N 2002bo Si appears to be mixed in 
with O (Stehle et al. 2005). The strong stratification observed 
in SNR 0519-69.0 suggests, therefore, that its supernova did not 
belong to the HVG class. 

This does not mean, however, that Kepler or Tycho did be- 
long to the HVG class, because other mechanisms may have 
destroyed the strong layering in the SNR phase. Kepler is in 
that sense a very interesting case, as the northern and south- 
ern regions show different behaviors. In the north the blast wave 
of Kepler appears to be interacting with a pre-existing shell of 
about 1 M , which caused a strong deceleration of the blast wave 
( Vink 2008 ). The lack of stratification in that region may there- 
fore be attributed to the strong deceleration, squeezing the dif- 
ferent layers closer together. 

Interestingly, the expansion rate of Tycho in the eastern part, 
where the ejecta are very close to the shock front, is also much 
lower than expected for a SNR expanding in a homogeneous 
medium (Kats uda et al.|2010| ). It is, therefore, possible that also 
in Tycho the initial stratification is no longer visible, as the lay- 
ers have been squeezed together due to the presence of a strong 
density gradient, and perhaps even mixed due hydrodynamical 
instabilities. 

Another explanation for the proximity of the ejecta to the 
shock front in Tycho is the effect of efficient cosmic ray accel- 
eration (Warren et a l.|2005| ). The efficient acceleration alters the 
equation of state due to the large pressure fraction attributed to 
relativistic particles and energy losses due to escape of the high- 
est energy cosmic rays ( |Decourchelle et aT|2000 |Helder et al.| 
2009). The distance between contact discontinuity and shock 
front is smaller in the east, whereas the X-ray synchrotron rims, 
a signature of cosmic ray acceleration, are more pronounced in 
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the west. Moreover, one should not confuse the distance between 
contact discontinuity and shock front, with a lack of stratification 
within the ejecta. Only if the equation of state within the ejecta 
is changed, can one explain the lack of stratification by invoking 
cosmic rays. This could be the case, if also the reverse shock is 
efficiently accelerating particles. For Tycho there is no evidence 
that this is the case, but the reverse shock of the core collapse 
SNR Cas A has been proposed as a site of efficient cosmic ray 
acceleration ( |Helder & Vink|2008| ). 

In addition to cosmic ray acceleration it has been proposed 
that Rayleigh-Taylor instabilities have brought ejecta close to 
the shock front ([Chevalier et al.|1992]|Wang & Chevalier|2001 



Badenes et al.|2006| l. If this is indeed the case, it may also have 
destroyed the initial stratification of the ejecta. It is, however, 
not clear why this mechanism would have affected the stratifi- 
cation of Tycho but not of SNR 0519-69.0, or the southern part 
of Kepler. So perhaps it is a combination of factors that made 
Tycho and the north of Kepler less stratified, than 0519-69.0 and 
the south of Kepler. 

5.3. The chemical composition of SNR 0519-69.0 

From our spectral analysis we derive ejecta masses for O, Si/S, 
Ar/Ca, and Fe. These masses should be regarded with some 
caution, because inhomogeneities in the density distribution can 
lead to some systematic errors. However, on face value, our es- 
timates indicate that the overall abundance pattern matches that 
for normal Type la supemovae, but that the oxygen mass is rather 
large, M * 0.2 - 0.3 M . 

According to a recent analysis of optical spectra of Type la 
supernovae ( |Mazzali et al. 2001) the amount of material from 
the C/O white dwarf burned into IME and more massive el- 
ements is remarkably homogeneous, corresponding to a total 
mass of roughly MrME.Fe.Ni ~ 1.1 M . This means that about 0.2- 
0.3 M consists of lower mass elements. This is consistent with 
the amount we infer for SNR 0519-69.0, provided that most of 
the outer layers consist of oxygen and that little unburned carbon 
is present. 

One of the differences between pure deflagration and delayed 
detonation explosion models is that the latter can burn material 
all the way to the outer edge of the white dwarf, leaving no car- 
bon. Optical spectroscopy of Type la supernovae reveal that little 
or no carbon is present, with the oxygen abundance a factor 100 
to 1000 times the carbon abundance (Mari on et al.| |2006). The 
relatively high mass of oxygen in 0519-69.0 is therefore consis- 
tent with optical spectroscopy of Type la supernovae, and sup- 
ports the idea that burning in the outer layers converts almost all 
carbon to oxygen. 

The amount of Fe that we infer, ~ 0.7 M Q , is also fairly typ- 
ical for Type la supernovae. Here we need to caution that not 
all the Fe may have been shock heated by the reverse shock, 
so the total iron mass may be higher. |Mazzali et al. ([2007) 
showed that the peak luminosity, for which they use the drop 
in magnitude within 15 days after maximum (Am 15) as a proxy, 
correlates with the total amount of iron/nickel synthesized (see 
also Woosley et al.|2007] >. Using their "Zorro diagram" we can 
translate our estimate of shocked iron (which includes decayed 
56 Ni), in an estimate for Amis- M Fe > 0.7 M Q corresponds to 
Am^ < 1.4. Combining the amount of IME that we estimate 
Mime = 0.2 - 0.4 M G with the "Zorro-diagram" suggests that 
Am 15 > 1.1. 

These constraints suggest that the supernova was a fairly nor- 
mal Type la event, not as bright as SN1991T (Am l5 = 0.95), 
but similar to for example the template Type la SN 1996X 



(Am^ = 1.3, Salvo et al. 2001). This is unlike the youngest 
known Type la SNR in the LMC, 0509-67.5, which was accord- 
ing to the light echo spectrum of the supernova an exceptionally 
bright Type la ( |Rest et al.|2008] ), similar to SN 1991T. 

There is some evidence that bright Type la supernovae are 



associated with younger stellar populations (Gallagher et al. 
|2008| ), but the starforming history around SNR 050 9-67.5 is 
not significantly younger than around SNR 0519-69.0 (Badenes 
|et al.|[2009| . This study, therefore, does not provide evidence 
for a one to one relation between population age and supernova 
brightness. In addition, |Howell et al. ([2009 ) show that the rela- 
tion between the mean galaxy stellar population age and Type 
la brightness has significant scatter. This may be no suprise, if 
the brightness of a Type la supernova is related to the properties 
and history of the progenitor white dwarf, whereas the time de- 
lay between birth of the binary system and the supernova event 



is also gov erned by the evolution of the secondary star ( Howell 
et^D^OO^n 



5.4. Analytical model 

From the measured locations of the reverse and forward shocks, 



we can derive typical velocities of the plasma, using the Truelove 
|& McKee] ( |1999| ) analytical models, with the ejecta density 
power n = 7. Setting the density of the circumstellar matter to 
"csm = 2 cirT 3 , supernova explosion energy to E — 10 51 erg, 
and the mass of the progenitor to Mwd = 1-4 M Q , we can plot the 
forward to reverse shock radii ratio versus the time that passed 
since the explosion. We see from the left panel of the Fig. 10 that 
these parameters correspond to age of the remnant of 460 ±150 



years, which is somewhat lower than the value estimated by Rest 



et al. (2005 1. The leftmost bottom panel of Fig. 10 shows that in 



this model, the ratio of the reverse to forward shock velocity is 
vrs/vbw - 1-2. 

5.5. Numerical HD models 



We run several hydrodynamical models (Soro kina et al.| |2004 
Kosenko 2006 1, in order to reproduce the observed features 
of SNR 0519-69.0. We set the CSM density to p C SM = 
3 x 10~ 24 gcm~ 3 . We considered two "toy" explosion mod- 
els by [Woosley et aL ( 2007 1 of thermonuclear supernovae: 
a "mildly mixed" c050403 (M Ni = 0.5M G , M Fe = 
0.4M G , M mE = 0.3M G ) and a "moderately mixed" m060403 
(M Ni = 0.6M G , Mp e = 0.4M G , M IME = 0.3M G ). The distribu- 
tion of the elements in these models is plotted in the right lower 
panels of Fig. [10] Velocity and density profiles of the models at 
the age of 560 years, which are based on these two initial set ups, 



are presented in Fig. 10 (top graphs) 



In the hydrodynamical models the contact discontinuity is 
further away from the blast wave and closer to the reverse shock 
as compared to the observed locations in the SNR 0519-69.0. 
Suggesting that some of the effects that may have brought the 
ejecta in Tycho's SNR close to the forward shock, whether 
pure hydrodynamical effects, or efficient cosmic ray accelera- 
tion, may to a lesser extent also have affected SNR 0519-69.0. In 
the near future we intent to produce models that incorporate cos- 
mic ray acceleration into the hydrodynamic models (Kosenko et 
al. in preparation). Note that Cassam-Chenai et al. (2008 ) found 



1 Nevertheless, some correlation between time delay and white dwarf 
properties can be expected, since a long time scale for the evolution 
of the primary star (the white dwarf that will explode), automatically 
means that the secondary star evolves on an even longer time scale. 
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Fig. 10. Left two plots: [Truelove & McKee ( 1999 1 analytical solution for SNR 0519-69.0. Radii ratio vs. the age (top panel); red 
dash-dotted lines show the measured value from the Chandra measurements (Sect. 3.1 1. Reverse shock velocity (in the frame of the 
ejecta) to forward shock velocity ratio vs. the age (bottom panel); red dash-dotted lines show values of age derived from the top 
plot. Right four plots: density, velocity (top panels) and abundances (bottom panels) profiles of the numerical models. In this set, 
left plots show the "mildly mixed" c050403 model, right plots — "moderately mixed" m060403. 



that for SN 1006 the contact discontinuity is also close to the 
shock front, to a lesser extent even for those regions that do not 
emit X-ray synchrotron radiation, where efficient cosmic ray ac- 
celeration is likely not important. Their models suggest that cos- 
mic ray acceleration alone is not sufficient to explain the obser- 
vations. This suggests that the distance of the contact continuity 
to the shock front is the result of a mixture of efficient cosmic ray 
acceleration and hydrodynamical instabilities (see also |Miceli| 
|et al.|2009] and the discussion in Sect. |5.2) >. 



and not an O/Ne/Mg white dwarf, which is consistent with type 
la explosion models (e.g. Woosle y et al.|2007) . 



5.6. RGS line broadening and the shock velocity 

Using the measured value of the emission lines broadening 
(cr RGS ), we can estimate the velocity of the forward shock. The 
temperature of each specie behind the reverse shock is calculated 
from (no temperature equilibration between species assumed) 



Estimated from our numerical models, the amount of swept 
up supernova matter is as follows: Mo - 7%, Ms;/s = 
12% 11%, M Ar /Ca = l%/2%, M Fe = 63% for c050403 and 
M = 3%, M si/S = 12%/7%, M Ar/Ca = l%/2%, M Fe = 71% 
for m060403 explosion model. With a total mass of the shocked 
ejecta of 1.2 M Q in both cases. We see that the fractional mass 
of Fe in c050403 model matches reasonably well the observed 
fractional mass as determined in Sect. |4.3| but the observed frac- 
tional mass of oxygen appears to be higher than indicated by the 
models. This may hint that the progenitor was an oxygen-rich 
white dwarf. The lack of neon in the RGS spectra and magne- 
sium in the EPIC spectra still argues in favor of C/O white dwarf, 



16 



(3) 



here m, — mass of species i, vrs — reverse shock velocity. Thus, 
the Maxwellian thermal broadening of the lines is 



0"thej 




V3 



-Vrs 



(4) 



As it is shown above, according to the analytical |Truelove &| 
|McKee| (19991 solution, for the case of SNR 0519-69.0 the rela- 
tion vrs = 1.2 vps holds. Also, according to the numerical mod- 
els, the value of the bulk velocity in the ejecta is by ~ 50% 
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lower than the bulk velocity in the shocked CSM (the velocity 
profiles of the HD models on the top right plots of Fig. 10 1, i.e. 
Obuikej - 0.5cr bulkcSM = 0.5 (3/4) vfs- Thus 



0"RGS 



thej 



+ CT 



bulk ej 



0.64 vps 



(5) 



Substituting the measured value of ctrgs = 1873 + 50 kms -1 
into eq. (|5jl, we obtain vps = 2927 ± 500 km s _1 . Where the error 
includes an estimate of the systematic uncertainties, which come 
from relation of ctrgs to vps- 

From the optical observation of the remnant we have inde- 
pendent measurements of the blast wave velocity, but these mea- 
surements give ambiguous values, [van Ade lsberg et al. (2008) 
modeled hydrogen spectra of the SNR and found shock veloci- 
ties of v FS = 2984+™ km s" 1 and v FS = 1178^ kms" 1 based 
on the observations of |Ghavamian et aL] ( |2007[ ) and |Smith et al.| 
(|1991|l respectively . Measurements of Ha-line along the rim by 
Tuohyetal.|(|l982| provide the value of v FS = 2800±300 km s _1 . 



Our estimate is the closest to the estimates by |Tuohy et al.| ( |1982) 
and |Ghavamia n et al.| ([2007 ). The lower value of the forward 
shock velocity in the measurements by Smith et al. ( 1991 1 is 



probably caused by the slit location they used, which intersects 
the bright eastern rim of the SNR, where the blast wave proba- 
bly encounters a dense cloud of interstellar matter and may have 
slowed down (Ghavamia n et al.|2007) . 

6. Conclusions 

We analyzed X-ray data of SNR 0519-69.0, which were ob- 
tained by both the XMM-Newton and Chandra observatories. 
We treated the spectra with spex fitting software and applied non- 
equiliblium ionization (NEI) models. Fitting of the data with the 
single-ionization timescale NEI model allowed us to measure 
average electron temperature of 3 keV, ionization timescale of 
2.3 x 10 10 s cm" 3 . 

The stratification of the elements is particularly pronounced 
in the radial emissivity profiles of the Chandra images of the 
SNR. Chandra ACIS spectra of specific regions of the SNR also 
points to the separation of the elements in the shell. This in- 
spired us to use a multi-component fitting model for the EPIC 
MOS and RGS spectral analysis, in which the observed layers 
are each represented by an NEI component dominated by one 
or two elements. This multi-component model provides an ade- 
quate description of both the broad band EPIC MOS spectra and 
the high spectral resolution RGS spectra in the 0.5-1.1 keV band. 
In addition it allowed us to characterize the Chandra ACIS spec- 
tra of several distinct regions. Using this approach we estimate 
the density of the circumstellar matter of 2.4 ± 0.2 cm -3 . 

The results of our analysis show that the X-ray data of 
SNR 0519-69.0 are consistent with mildly or moderately mixed 
delayed detonation explosion models, which result in clear sep- 
aration of the outermost oxygen layer (produced during the fi- 
nal detonation phase), followed inwards by the partially mixed 
inner layers of intermediate mass elements and iron (deflagra- 
tion phase of the explosion). The mass fraction that we infer are 
M * 32%, M si/S ~ 7%/5%, M Ar+Ca ~ 1% and M Fe * 55%, 
which should be considered with some caution, as we had to 
infer a volume filling factor to convert from emission measure 
to mass. Taken at face value, these mass fraction are consistent 
with recent observational results of Type la supernovae (Mazzali 
|et al.|2007] ), which show that all material out to an enclosed mass 
of ~ 1.1 Mq is burned into IME and heavier elements, leaving 
about 0.2 - 0.3 M Q for oxygen, provided all carbon was burned. 



The mass fractions of IME and Fe in SNR 0519-69.0 puts con- 
traints on the sub-classification of the supernova, suggesting a 
Type la event of normal peak luminosity. 

The RGS data allowed us to measure the line broadening of 
the shocked ejecta of 1873 + 50 kms~ 3 , which helped to con- 
strain the dynamical properties of the remnant and estimate the 
forward shock velocity of 2770 ± 500 km s _1 . This value is con- 



sistent with the measurements reported by Tuohy et al. ( 1982 1 
and Ghavamia n et al.| ( |2007| l, but inconsistent with |Smith et aL] 

qi99i) . 

Acknowledgements. We thank the anonymous referee for his/her insightful 
comments, which helped to broaden the scope of this paper. We also thank 
S.I. Blinnikov for providing us with the thermonuclear explosion models and 
also for the fruitful discussions on explosion mechanisms. DK is supported by 
an "open competition" grant and EH and JV are supported by a Vidi grant from 
the Netherlands Organization for Scientific Research (PI J. Vink). The work is 
also partially supported in Russia by RFBR under grant 07-02-00961 . 



References 

Anders, E. & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 
Badenes, C, Borkowski, K. J., Hughes, J. P., Hwang, U., & Bravo, E. 2006, 

Astrophys. J., 645, 1373 
Badenes, C, Bravo, E., Borkowski, K., & Inmaculada, D. 2003, Astrophys. J., 

593,358 

Badenes, C, Hams, J., Zaritsky, D., & Prieto, J. L. 2009, Astrophys. J., 700, 727 
Benetti, S., Cappellaro, E., Mazzali, P. A., et al. 2005, Astrophys. J., 623, 1011 
Blinnikov, S. I., Ropke, F. K., Sorokina, E. I., et al. 2006, A&A, 453, 229 
Branch, D., Baron, E., Hall, N., Melakayil, M, & Parrent, J. 2005, PASP, 117, 
545 

Branch, D., Dang, L. C, & Baron, E. 2009, PASP, 121, 238 
Cassam-Chenai, G., Decourchelle, A., Ballet, J., et al. 2004, A&A, 414, 545 
Cassam-Chenai, G, Hughes, J. P., Reynoso, E. M., Badenes, C, & Moffett, D. 

2008, Astrophys. J., 680, 1 180 
Chevalier, R. A., Blondin, J. M., & Emmering, R. T. 1992, Astrophys. J., 392, 

118 

Decourchelle, A., Ellison, D., & Ballet, J. 2000, Astrophys. J., 543, 57 

den Herder, J. W., Brinkman, A. C, Kahn, S. M., Branduardi-Raymont, G., & 

Thomsen, K. e. a. 2001, A&A, 365, 7 
Gallagher, J. S., Garnavich, P. M., Caldwell, N., et al. 2008, Astrophys. J., 685, 

752 

Ghavamian, P., Blair, W. P., Sankrit, R., Raymond, J. C, & Hughes, J. P. 2007, 

Astrophys. J., 664, 304 
Hamilton, A. J. S., Fesen, R. A., & Blair, W. P. 2007, MNRAS, 381, 771 
Helder, E. A. & Vink, J. 2008, Astrophys. J., 686, 1094 
Helder, E. A., Vink, J., Bassa, C. G., et al. 2009, Science, 325, 719 
Howell, D. A., Sullivan, M., Brown, E. F, et al. 2009, Astrophys. J., 691, 661 
Hughes, J. P., Hayashi, I., Helfand, D., et al. 1995, ApJ, 444, 81 
Kaastra, J., Mewe, R., & Nieuwenhuijzen, H. 1996, in UV and X-ray spec- 
troscopy of astrophysical and laboratory plasmas (K. Yamashita and T. 
Watanabe), 411 

Katsuda, S., Petre, R., Hughes, J. P., et al. 2010, Astrophys. J., 709, 1387 
Koerwer, J. F. 2009, AJ, 138, 1 
Kosenko, D. 2006, MNRAS, 369, 1407 

Kosenko, D., Vink, J., Blinnikov, S., & Rasmussen, A. 2008, A&A, 490, 223 
Lewis, K. T., Burrows, D. N., Hughes, J. P., et al. 2003, Astrophys. J., 582, 770 
Lucy, L. B. 1974, AJ, 79, 745 

Marion, G. H., Hoflich, P., Wheeler, J. C, et al. 2006, Astrophys. J., 645, 1392 
Mazzali, P. A., Ropke, F. K, Benetti, S., & Hillebrandt, W. 2007, Science, 315, 
825 

Miceli, M., Bocchino, F, Iakubovskyi, D., et al. 2009, A&A, 501, 239 
Panagia, N. 2005, II Nuovo Cimento B, 120, 667 

Perlmutter, S., Aldering, G, Goldhaber, G., et al. 1999, Astrophys. J., 517, 565 
Phillips, M. M., Lira, P., Suntzeff, N. B., et al. 1999, AJ, 1 18, 1766 
Pietrzynski, G, Thompson, I. B., Graczyk, D., et al. 2009, Astrophys. J., 697, 
862 

Rest, A., Matheson, T, Blondin, S., et al. 2008, Astrophys. J., 680, 1 137 
Rest, A., Suntzeff, N. B., Olsen, K, et al. 2005, Nature, 438, 1 132 
Richardson, W. H. 1972, Journal of the Optical Society of America, 62, 55 
Riess, A. G., Filippenko, A. V, Challis, P., et al. 1998, AJ, 1 16, 1009 
Ropke, F. K. & Bruckschen, R. 2008, New Journal of Physics, 10, 125009 
Russell, S. C. & Dopita, M. A. 1992, Astrophys. J., 384, 508 
Salvo, M. E., Cappellaro, E., Mazzali, P. A., et al. 2001, MNRAS, 321, 254 



12 



Kosenko et al.: Kinematics and chemical stratification of SNR 0519-69.0 



Smith, R. C, Kirshner, R. R, Blair, W. P., & Winkler, R R 1991, Astrophys. J., 
375, 652 

Sorokina, E., Blinnikov, S., Kosenko, D., & Lundqvist, R 2004, Astronomy 
Letters, 30, 737 

Stehle, M., Mazzali, R A., Benetti, S., & Hillebrandt, W. 2005, MNRAS, 360, 
1231 

Truelove, J. K. & McKee, C. F. 1999, ApJS, 120, 299 

Tuohy, I. R., Dopita, M. A., Mathewson, D. S., Long, K. S., & Helfand, D. J. 

1982, Astrophys. J., 261, 473 
Turner, M. J. L., Abbey, A., Arnaud, M., Balasini, M., & Barbera, M. e. a. 2001, 

A&A, 365, 27 

van Adelsberg, M., Heng, K., McCray, R., & Raymond, J. C. 2008, Astrophys. 

J., 689, 1089 
Vink, J. 2008, Astrophys. J., 689, 231 

Vink, J., Laming, J. M., Gu, M. R, Rasmussen, A., & Kaastra, J. S. 2003, 

Astrophys. J., 587, L31 
Wang, C.-Y. & Chevalier, R. A. 2001, Astrophys. J., 549, 1 1 19 
Warren, J. S., Hughes, J. P., Badenes, C, et al. 2005, Astrophys. J., 634, 376 
Williams, R., Petre, R., & Holt, S. S. 2001, AIP Conference Proceedings, 565, 

185 

Willingale, R., West, R. G., Pye, J. P., & Stewart, G. C. 1996, MNRAS, 278, 749 
Woosley, S. E., Kasen, D., Blinnikov, S., & Sorokina, E. 2007, Astrophys. J., 
662,487 



